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The  breaking  cf  the  H-H  bond  in  H3  and  a  C-H  bond  in  CH4  on 
both  discrete  transition  metal  complexes  and  on  Ni  and  Ti  surfaces  is  studied, 
and  the  essential  continuity  and  similarity  of  the  physical  and  chemical  processes 
in  the  two  cases  is  demonstrated.  We  begin  with  an  orbital  analysis  of  oxida¬ 
tive  addition,  delineating  four  basic  interactions:  H-H  or  C-H  c  -»  M  elec¬ 
tron  transfer,  the  reverse  M  -»  a*  transfer  (both  weakening  the  c  bond,  form¬ 
ing  the  M-H  bond) ,  a  repulsive  interaction  between  o  and  metal  filled  orbitals , 
and  a  rearrangement  of  electron  density  at  the  metal.  The  molecular  cases  anal¬ 
yzed  in  detail  are  d®  MLs  ,  d8ML4  and  CpM'L.  Coordinative  unsaturation  is 
necessary,  and  consequently  a  •*  M  electron  transfer  dominates  the  early 
stages  of  the  reaction.  Steric  effects  are  important  for  CH4  reaction.  Ac¬ 
tivation  in  angular  MI*  or  CpM'L  is  achieved  through  a  destabilized  yz 
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MO,  and  d10ML3,  MLS  candidates  for  activation  are  described.  For  our 
study  of  the  surface  we  develop  tools  such  as  projections  of  the  density  of 
states  and  crystal  orbital  overlap  populations  -  the  extended  structure  ana¬ 
logues  of  a  population  analysis.  These  allow  a  clear  understanding  of  what 
happens  when  an  Ha  or  a  CH*  molecule  approaches  a  surface.  Because  of 
the  higher  energy  of  the  occupied  metal  orbitals  on  the  surface  the  M  -+  o* 
interaction  leads  the  reaction.  There  are  great  similarities  and  some  differ¬ 
ences  between  the  activation  acts  in  a  discrete  complex  and  on  a  surface. 
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In  this  paper  we  will  try  to  understand  how  an  H-H  or  C-H  bond 
can  interact  and  eventually  break  in  the  proximity  of  one  or  more  transition 
metal  centers.  We  analyze  this  problem  both  for  discrete  complexes  and  for 
a  clean  metal  surface;  indeed  the  most  interesting  aspect  of  our  study  will 
be  the  comparison  of  similarities  and  differences  between  the  chemistry  that 
goes  on  in  an  inorganic  complex  and  on  a  metal  surface. 

Let  us  review  the  experimental  background  of  this  problem.  Until  re¬ 
cently  there  was  a  nice  sharp  dichotomy  in  the  chemistry  of  Ha  with  trans¬ 
ition  metal  complexes .  If  Ha  interacted  at  all,  it  reacted  completely,  yield¬ 
ing  in  an  oxidative  addition  process  a  metal  dihydride,  l.  This  species  was 
sometimes  observed,  more  often  inferred  as  it  was  consumed  rapidly  in  some 
subsequent  rapid  chemistry.  *  Recently  the  first  well-characterized  Ha  com- 
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plex  was  observed.  This  is  Z>  a  side-one  bonded  complex  with  a  d  ML5 

fragment.  The  H-H  distance,  available  from  a  neutron  diffraction  study,  is 
0.75  ±  0.16  A. 
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For  the  interaction  of  a  C-H  bond  with  one  or  more  metal  atoms  the 
experimental  history  is  much  richer.  Over  the  years  it  has  become  apparent 
that  a  C-H  a  bond  can  interact  in  a  bonding  way  with  a  coord inatively  un¬ 
saturated  metal  center  (16  or  less  electrons  around  the  metal)  and  in  so 
doing  allow  the  metal  to  achieve  or  approach  the  stable  18  electron  configu- 
ration.  Unsaturation  at  the  metal  and  proximity  are  required.  The  inter- 

molecular  cases,  most  of  which  are  quite  recent,  proceed  on  to  oxidative 
4 

addition,  3. 
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The  intramolecular  examples,  ones  in  which  the  interacting  alkyl 
group  is  somehow  tethered  to  the  metal  atom,  have  been  revealing  in  show¬ 
ing  us  details  of  the  initial  states  of  metal-CH  interaction.  There  is  by 
now  an  ample  store  of  structural  or  spectroscopic  evidence  for  intramolecular 
M-CH  interaction  with  a  variety  of  geometries,  coordination  numbers  and 
electron  counts  at  one  or  more  metal  atoms.  Precise  structure  determinations, 
utilizing  neutron  diffraction ,  show  short  M  to  H  (and  C)  contacts,  and  un¬ 
usually  long  C-H  bonds  (the  world's  record  now  stands  at  1.19A^).  C-H 
stretching  frequencies  often  are  dramatically  lowered  and  C-H  coupling 
constants  as  well.  There  is  evidence  of  both  linear,  4a  and  triangular,  4b. 
interaction  geometries. 
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Hs  and  alkanes  are,  of  course,  chemisorbed, dissociated, and  reassembled 
on  many  transition  metal  surfaces  of  varying  degrees  of  cleanliness.  In  recent 
times  the  reactivity  of  definite  crystal  planes  has  been  studied  in  some  detail, 
and  we  are  beginning  to  gain  information  on  the  microscopic  structure  of  the 
product  surface.** 

Theoretical  studies  of  surfaces  and  their  interaction  with  molecules  are 

now  being  done  by  several  groups.  We  want  to  single  out  for  special  mention 

7 

here  the  work  of  R.  Baetzold,  E.  Shustorovich  and  E.  Muetterties  because  it 
anticipates  many  of  the  conclusions  that  we  eventually  reach  about  the  surface 
in  particular  concerning  the  role  of  the  substrate  o*  orbitals  and  the  direction 
of  electron  flow  during  surface- substrate  interaction.  Other  theoretical  studies 
will  be  mentioned  in  the  course  of  the  paper. 
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Charge  Transfer  and  Bond  Making  in  Oxidative  Addition 

At  the  risk  of  repeating  what  is  obvious  let  us  examine  the  essential 

features  of  oxidative  addition,  correlating  the  basic  ideas  of  electron  transfer 

and  oxidation,  reduction  with  the  way  in  which  these  appear  in  a  molecular 

8 

orbital  description  of  a  process. 

Drawing  5  is  a  schematic  illustration  of  the  level  transformation  in  a 


transition  metal  complex  reacting  with 


- Acc 

■H —  Don 


Hs  .  The  MLn  complex  is  represented  \ 
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by  a  band  of  occupied  levels  and  a  band  of  unoccupied  ones.  The  metal 
surface  will  be  no  different.  One  of  the  metal  filled  levels  (Don  for  Donor) 
and  one  of  the  unfilled  ones  (Acc  for  Acceptor)  is  singled  out,  for  reasons 
which  will  soon  become  apparent.  At  the  end  of  the  reaction  two  new  M-H  o 
bonds  form,  and  of  course  their  corresponding  antibonding  combinations.  In 
the  conventional  Wernerian  scheme  of  counting  ligands  as  two-electron  o 
donors  the  four  electrons  of  the  two  new  M-H  c  bonds  are  assigned,  for 
electron  counting  purposes,  to  the  ligands,  H  .  It  is  this  convention  which 
makes  the  metal  go  from  a  dn  to  a  dn  2  electron  count,  and  makes  us  call 
this  reaction  an  oxidative  (at  the  metal)  addition. 

Formalisms  are  convenient  fictions  which  contain  a  piece  of  the  truth  -  and 
it  is  so  sad  that  people  spend  a  lot  of  time  arguing  about  the  deductions  they  draw, 
often  ingeniously  and  artfully,  from  formalisms,  without  worrying  about  their 
uncerlying  assumptions.  The  "complex"  or  dative  bonding  picture  which  led  to 
"oxidation  at  metal"  of  course  is  an  exaggeration.  The  M-H  c  bonds  are  in 
good  pari  covalent.  To  the  extent  that  they  are  so,  the  real  d  electron  population 
at  the  metal  moves  back  from  d°  ^  toward  dn.  To  the  extent  that  it  probably 
never  quite  gets  back  to  dn  it  is  still  informative  to  call  this  an  oxidative  addition. 

What  the  "oxidative  addition"  formalism  conceals  and  a  molecular 
orbital  picture  reveals  is  that  in  the  course  of  this  reaction  there  has  to  be 
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a  two  way  flow  of  electron  density,  from  the  metal  to  the  new  ligands  and 
in  the  reverse  direction. 

Consider  the  M-H  o  bonds  in  the  product.  In  a  localized  represen¬ 
tation  they  are  shown  in  6,  and  the  equivalent  delocalized  picture  in  7. 

6 

The  delocalized  orbitals  are  labelled  as  S  or  A  according  to  their  symmetry 
or  antisymmetry  with  respect  to  the  two-fold  axis  or  mirror  that  interchanges 
them. 

Where  did  7S  and  7A  originate?  They  came  from  the  interaction  of 
aHa  with  metal  Acc  and  metal  Don  with  o*Hg  .  This  is  shown  in  £  and  in 


another  way,  focussing  on  the  evolution  of  the  orbitals ,  in  9.  What  g  shows 
clearly  is  the  two  way  charge  transfer  and  the  coupling  of  electron  transfer  and 
bonding  changes .  The  symmetric  M-H  combination  evolves  from  cHa  by 
mixing  in  of  a  metal  acceptor  orbital  of  appropriate  symmetry.  The  result  is 
electron  transfer  from  Hac,  decreased  H-H  bonding  and  increased  M-H 
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bonding.  The  A  combination  electron  transfer  is  in  the  opposite  direction, 
for  this  orbital  is  originally  on  the  metal.  Electron  transfer  to.  Ha^*  has 
as  a  consequence  decreased  H-H  bonding  and  increased  M-H  bonding. 

Note  that  both  these  interactions  lead  to  H-H  bond  weakening  and 
M-H  bond  formation,  even  though  they  accomplish  these  actions  by  charge 
transfer  in  different  directions. 

The  molecular  orbital  description  makes  it  clear  that  when  oxidative 
addition  is  complete  there  must  have  occurred  electron  transfer  from  metal  to 
Hs  or  RH  and  in  the  reverse  direction.  But  there  is  no  requirement  that 
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the  electron  flow  be  balanced  at  every  stage  of  the  reaction.  In  fact  the 
experimental  evidence  for  the  requirement  of  coordinative  unsaturation  of 
the  metal  in  activation  on  discrete  complexes  makes  it  clear  that  in  these 
molecules  the  important  initial  electron  flow  is  from  Ha  to  metal.  As  we 
will  see  metal  surfaces  may  be  different. 

So  far  we  have  identified  the  two  most  important  bonding  interactions 
between  hydrogen  molecule  or  an  alkane  and  a  discrete  transition  metal 
complex  or  a  metal  surface.  They  are  repeated  in  JJD,  labeled  as  and 
^2 ^  now.  These  are  two-orbital  two-electron  bonding  interactions.  Two 
further  interactions  must  be  thought  about.  Interaction  (^J  is  the  two-or¬ 
bital  four-electron  perforce  destabilizing  interaction  between  filled  orbitals 
of  substrate  and  surface  (or  complex).  It  is  in  this  interaction  that  one-elec¬ 
tron  theories  of  the  extended  Hiickel  type  find  what  chemists  normally  call 


"steric  effects" .  Interaction 


© 


is  destabilizing,  and  leads  to  some  M-H 


antibonding.  It  is  the  primary  source  of  barriers  to  C-H  activation. 
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In  addition  to  these  three  interactions, which  operate  for  an  MLn  com¬ 


plex  as  well  as  for  a  metal  surface,  there  is  another  interaction 


■0 


which 


is  generally  important  only  for  metal  surfaces,  where  there  are  many  closely 
spaced  levels.  What  happens  on  the  surface,  as  we  will  see  in  great  detail 
below,  is  that  some  levels,  more  localized  on  the  surface  than  in  the  bulk, 
and  even  on  the  surface  distinguished  by  reaching  out  toward  the  substrate, 
some  levels  interact  to  a  greater  extent  than  others.  They  do  so,  of  course, 
through  the  primary  interactions  3©  «ndQ.  But  since  in  a  solid  or 
a  surface  levels  are  closely  spaced  around  the  Fermi  level,  the  net  result 
of  such  primary  interactions  of  substrate  and  "surface  states"  is  a  shift  of 


electron  density  between  bulk  and  surface,  and  even  within  the  surface.  This 


interaction  is  poorly  represented  in  H)  by 


0 


but  its  significance  will 


eventually  become  clearer  as  we  describe  it  in  more  detail  later. 


We  are  now  ready  to  proceed  with  an  analysis  of  several  specific  cases, 

to  see  these  interactions  in  action.  But  first  let  us  describe  the  computational 

9 

methodology  we  use.  This  is  the  extended  Huckel  method,  with  particulars 
described  in  Appendix  1.  There  is  a  special  problem  which  this  transparent 


and  simple  procedure  brings  with  it.  The  method  is  not  reliable  for  bond 
distance  changes,  and  Ha  in  particular  is  a  pathological  case  in  which  the 
two  atoms  collapse.  So  the  study  of  potential  energy  surfaces  where  H-H 
or  C-H  bonds  are  made  or  broken  would  seem  to  be  an  inappropriate  applica- 


i 


tion  of  the  extended  Huckel  method.  In  fact  this  is  so,  and  since  we  cannot 


trust  the  method  for  bond  distances  we  do  not  calculate  complete  potential 
energy  surfaces.  Instead  we  limit  ourselves  in  general  to  the  study  of  select  ap¬ 
proaches  -  for  instance  an  Ha  coming  on  parallel  or  perpendicular  to  a  surface  - 
and  focus  on  that  aspect  of  the  electronic  redistribution  which  the  extended  Huckel 
method  from  our  experience  is  likely  to  get  right.  This  is  the  magnitude  and 
nodal  character  of  orbital  interactions. 

We  also  apply  consistently  the  language  and  formalism  of  simple  per¬ 
turbation  theory,  in  particular  the  second  order  expression  for  the  interaction 
of  two  levels: 
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Extended  Huckel  arguments,  especially  in  the  fragment  orbital  analysis, 
translate  directly  into  perturbation  arguments.  It  is  this  combination  of  ex¬ 
tended  Huckel  calculations,  and  perturbation  theory  based  thinking  within  a 
one-electron  frontier  orbital  picture  that  makes  us  feel  more  sanguine  about 


the  results  of  what  would  otherwise  have  been  a  pretty  unreliable  calculation. 


Ha  and  CH*  as  substrates 


The  orbitals  of  both  molecules  are  familiar.  Within  a  simple  single 
configuration  picture  the  valence  orbitals  are  filled  Cg  in  Ha  ,  an  ai  +  ta 
set  in  CH*  ,  and  the  corresponding  unfilled  cu*  in  Ha  and  ai  *  and  t3  * 
in  CH4  .  The  orbital  energies  as  given  by  the  extended  Huckel  method  are 
shown  in  Figure  1.  The  C-H  bonding  in  CH*  is  distributed  over  the  ai  and 
t3  set,  but  is  mainly  in  the  ta  component.  If  we  focus  on  that  orbital  as  the 
C-H  bond  and  then  compare  CH4  and  Ha  then,  as  far  as  energetics  are  con¬ 
cerned,  the  two  molecules  are  equally  good  (poor)  acceptors,  but  CH*  is  a 
better  (but  still  not  good)  donor,  as  its  ta  set  is  some  two  eV  higher  in 
energy  than  Ka  a. 

Figure  1  here 

The  numerator  of  the  perturbation  sum,  i  H^j  i3,  is  not  to  be  forgotten. 
Coefficients  of  the  relevant  orbitals  are  given  in  H.*  Note  first  the  spectacular 
difference  between  the  ag  and  cu*  H  coefficients  in  Ha  .  This  is  a  result 
of  including  the  overlap  in  the  normalization  of  the  molecular  orbitals .  An 
immediate  consequence  is  that  c*  orbitals,  acting  through  the  numerator 
of  the  perturbation  expression,  will  have  much  more  "power"  in  the  in¬ 
teraction.  This  will  compensate  sometimes  (as  we  will  see,  especially 


1 


Energy  (eV) 


-  15a  - 


-  16  - 


on  transition  metal  surfaces)  for  their  very  high  energy,  which,  acting 
through  the  denominator  of  the  perturbation  expression,  makes  it  difficult 
for  them  to  have  much  influence.  Note  the  similarity  of  the  effective  H 
coefficients  in  CH4  and  Ha  . 

We  are  now  ready  to  proceed  with  calculations  on  the  addition  of  these 
molecules  to  various  MLn  fragments.  The  reader  should  note  that  our  calculations 
are  not  the  only  ones  extant,  and  that  several  others  have  been  published. 
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A  Prototype  Mononuclear  Transition  Metal  Fragment,  Cr(CO)6 

Why  Cr(CO)s  ?  Complexes  of  Cr(CO)s  with  methane  and  hexane  in 
in  low  temperature  matrices  have  been  detected.  ^  The  only  well-character¬ 
ized  Ha  complexes  M(CO)^PR3)aH2  ,  M=Mo,  W  are  closely  related  to 

Cr(CO)s2  ,  and  many  of  the  cited  intramolecular  cases  of  C-H  activation  can 

5X-t 

be  related  back  to  this  model. 

The  total  energy  of  an  Ha  frozen  at  H-H  0.74  A  approaching  a  C*v 
octahedral  fragment  Cr(CO)s  is  shown  in  Figure  2.  Two  approach  geometries 
were  studied,  a  "perpendicular"  and  "parallel"  mode.  These  are  sufficiently 

Figure  2  here 

common  in  the  subsequent  discussion  that  it  is  best  to  describe  them  more 
precisely  in  X2.  and  12.  •  "Perpendicular"  means  H-H  (or  eventually  C-H) 
colinear  with  the  metal  atom,  "parallel"  means  turned  by  90°  ,  so  that  both 
M-H  distances  are  equal.  The  LnM-Hs  separation,  somewhat  arbitrarily , 
is  defined  as  the  distance  to  near  hydrogen  in  the  perpendicular  geometry , 
but  to  the  H-H  centroid  in  the  parallel  form. 
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Fig.  2  Total  energy  along  the  perpendicular  and  parallel  approaches 
of  H£  to  Cr  (CO) g -  The  common  energy  zero  for  both  curves  is 
for  the  two  fragments  at  infinite  separation. 
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The  two  energy  curves  are  attractive  over  a  substantial  range  of  ap¬ 
proach  distances,  and  the  parallel  minimum  is  deeper.  This  is  consistent 
with  the  single  structure  known,  2,  but  why  is  it  so? 

A  level  interaction  diagram  (Figure  3)  is  illuminating.  Singled  out  in 
the  middle  of  the  construction  are  the  familiar  orbitals  of  a  16  electron  Cw 

Figure  3  here 


12 

MLs  fragment  -  a  t2g  set  below  a  well-directed  a1  hybrid.  In  the 
perpendicular  geometry  the  tag  set  is  untouched,  while  both  og  and  ou* 
interact  in  a  typical  three-orbital  pattern  with  ML5  ai  .  The  significant  bond¬ 
ing  mixing  is  of  type  in  10.»  between  og  and  a!  of  Cr(CO)s  in 
Figure  3.  It  gives  rise  to  stabilization,  but  the  stabilization  is  not  great  as 
it  is  underlain  by  a  repulsive  base  of  four-electron  destabilizing  interactions. 

In  the  parallel  geometry  the  interaction  between  og  and  ai  ,  though 
somewhat  different  in  spatial  configuration,  in  fact  is  not  significantly  worse 
in  overlap.  Now  a  strong  interaction  of  type  ^2),  metal  acting  as  a  donor 
(xz)  toward  the  ligand  as  acceptor  (au*),  is  allowed  by  symmetry,  whereas 
it  was  forbidden  in  the  perpendicular  approach.  There  is  electron  transfer 
from  a  (0.110  electrons  at  L=2.0  A)  and  into  c*  (0.032  electrons  at 
L=2 . 0  n) .  No  wonder  the  parallel  geometry  is  preferred. 

The  situation  is  changed  dramatically  if  the  H-H  bond  is  substituted 
by  C-H  of  CH4  (Figures  4  and  5).  The  perpendicular  approach  is  still 
attractive,  but  the  parallel  one  is  not,  becoming  strongly  repulsive.  The  level 


I 


Fig.  3  Interaction  diagram  for  Hg  and  Cr  (C0)g  for  a  perpendicular 

(left)  and  parallel  (right)  approach.  The  diagram  is  schematic 

In  the  position  of  the  a„  and  a*  levels  before  and  after  inter- 

g  u 

action. 
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analysis  reveals  that  the  origin  of  stabilization  in  the  perpendicular  approach 

Figures  4  and  5  here 

is  the  same  for  CH*  as  Ha  .  In  the  parallel  geometry  something  different 
is  happening, for  instead  of  xz  going  down  in  energy  (Ha),  in  the  case  of 
methane  both  xz  and  yz  go  up.  The  reason  for  this  behavior  may  be  seen 
from  the  fragment  overlaps  {at  La2.0  A)  in  14_. 


14  <xz  |  cr>  =  0.08 


The  Ha  xz  -  c*  overlap  is  big,  but  that  of  the  corresponding  o*  component 
of  CH4  is  small.  The  metal  d  orbital  is  mismatched  with  the  methane,  and 
samples  the  rear  of  the  CH  a*  combination.  Furthermore,  there  is  a  substantial 
overlap  between  metal  xz  and  the  occupied  C-H  o  orbital.  In  fact  this  re¬ 
pulsive  effect  dominates  pushing  xz  up  in  energy  (Figure  5)  as  it  interacts 
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Fig.  4  Total  energy  for  perpendicular  and  parallel  approaches  of 
CH^  to  Cr  ( CO ) 5 .  The  common  energy  zero  for  both  curves  is 
for  both  fragments  at  infinite  separation. 


1 


Fig.  5  Interaction  diagram,  schematic  at  the  top  and  bottom  of  the 
energy  range,  for  perpendicular  (left)  and  parallel  (right) 
approaches  of  CH^  to  Cr  (C0)g. 
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more  with  C-H  o  than  with  a*.  The  metal  yz  orbital  also  goes  up  as  a 
result  of  a  similar  four  electron  destabilization  with  another  member  of  the 
ta  set. 

What  we  have  is  the  dominance  of  two  orbital  four  electron  repulsions 
ever  the  attractive  bonding  forces.  We  think  that  we  can  identify  the  repulsive 
effects,  as  safely  as  one  can  do  it  within  the  framework  of  a  one  electron  theory, 
with  steric  effects . 

Thus,  when  a  methane  approaches  a  Cr(CO)5  in  a  parallel  geometry 
the  steric  problems  of  that  approach  overrule  the  favorable  geometric  arrange¬ 
ment.  Rotation  of  the  methyl  group  around  the  H-C  bond  does  not  alleviate 
the  trouble.  Can  a  geometry  intermediate  between  parallel  and  perpendicular 
achieve  a  compromise? 

We  tried  15.  and  J.6,  in  various  geometries .  Stabilization  was  achieved 


15 
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for  some  geometries,  for  instance  at  M-H  2.0  A,  .15  is  bound  for  all 
0  ^  130°  .  But  the  perpendicular  configuration,  0  =  180°  ,  is  most  stable. 
These  are  intermolecular  cases,  and  we  are  certain  that  in  special  intra¬ 
molecular  geometries ,  where  the  C-H  bond  is  so  suspended  near  the  metal 
that  steric  effects  are  minimized,  that  an  appropriate  triangular  geometry  with 
partial  M-H  and  M-C  bonding  is  attainable. 

A  particularly  interesting  example  of  intramolecular  interaction  occurs 
in  the  series  of  16  electron,  near  octahedral  complexes  of  type  17.  These 

L 


17 


exhibit  short  M’,#H  and  M**’C  contacts  associated  with  fluxional  behavior 
of  the  H  atom^a,t  V.  Previous  extended  Hiickel  calculations  on  the  model 
[Co(rf  -  alkenyl)(PH3)3]2+  have  pointed  out  a  low  lying  orbital  in  which  both 
M'^’H  and  M'^’C  interactions  are  bonding^.  Our  calculation  on  the 
isoelectronic  [Fe(,n3-  alkenyl)(CO)3 ]*  complex  shows  similar  results :  two 
relatively  delocalized  °CH  orbitals  of  the  methyl  group  are  stabilized  by  the 
metal  LUMO.  In  the  two  corresponding  bonding  MO's,  the  overlap  popu¬ 
lations  are  0.060  for  M-H  and  0.030  for  M-C.  The  corresponding  values 
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for  the  total  overlap  populations  are  0.075  and  0.  006  respectively .  Despite 
the  short  M-C  distance,  there  is  weak  interaction.  The  reason  is  that  the 
bonding  attraction  between  the  accepting  hybrid  orbital  and  the  bond  is 

balanced  by  a  secondary  H-C  antibonding  effect,  a  consequence  of  the  bond¬ 
ing  interaction  between  the  metal  orbital  and  the  polyenyl  it  system. 


l 


dJ  ML4  and  CpML  Systems  for  Activinq  C-H  Bonds 

Short  M-H  contacts  have  been  shown  to  exist  in  planar  or  near  planar 
16  electron  complexes  of  type  18^  ^ . 


18  19 


Although  the  accuracy  of  the  H  atom  positions  is  poor,  the  X-ray  crystallographic 
results  suggest  that  the  M-H  distances  lie  in  the  range  2. 3-3.0  a,3  _I'^  sub¬ 
stantially  longer  and  weaker  than  the  corresponding  distances  in  MLs  com¬ 
plexes.  More  interesting,  perhaps,  is  the  fact  that  we  can  place  in  the  d9ML4 
category  the  recent  exciting  cases  of  C-H  activation  using  the  CpML  fragment, 
M=Rh,  Ir,  19.4a~C 

Our  model  study  examined  the  approach  of  Hs  and  CH4  to  square  planar 


and  angular  Rh(C0)4*  fragments,  2j0  and  ZX. 
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Figure  6  shows  the  computed  total  energy  curves  for  Ha  .  Both  approaches  to 

the  square  planar  fragment  are  repulsive.  At  large  separations  the  perpendicular 

10a 

approach  is  preferred, in  agreement  with  previous  calculations  by  Sevin  and  by 

Figure  6  here 

10b 

Dedieu  and  Strich  .  When  the  square  planar  fragment  is  bent  back  to  the 

angular  one  ZX>  the  parallel  approach  is  greatly  stabilized. 

The  reasons  for  this  behavior  are  made  clear  by  examining  the  frontier 
13 

orbitals  of  £0  and  2X,  and  comparing  them  to  those  of  the  MLs  fragment, 
Cr(CO)s  .  This  is  done  in  ZZ •  The  crucial  difference  between  de  square 
planar  ML*  and  d  MLs  is  the  presence  of  the  occupied  z  in  the  former. 


'"'s 
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Fig.  6  Total  energy  for  perpendicular  and  parallel  approaches  of 
to  a  square  planar  and  angular  Rh  (CO)^.  The  common  energy 
zero  for  all  curves  is  at  infinite  separation  of  the  fragments. 
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and  the  very  different  makeup  of  the  LUMO.  The  za  orbital  introduces  an 
additional  high-lying  orbital  of  axial  symmetry,  capable  of  interaction  with 
a  and  a*  of  Ha  in  the  perpendicular  geometry,  with  o  in  the  parallel  one  . 

The  dominant  effect  is  the  four-electron  destabilizing  one,  and  this  is  what 
makes  the  two  square-planar  approaches  in  Figure  6  unfavorable.  The  acceptor 
orbital  of  ML*  is  a  poor  counterpart  of  the  MLe  hybrid.  The  ML*  orbital 
which  is  the  LUMO  is  only  15%  metal  p,  and  predominantly  ligand  ir*. 

Not  suprisingly,  given  what  we  know  about  the  chemistry  of  square-planar  d8 
complexes ,  the  coordinative  unsaturation  or  acceptor  power  of  square-planar 
ML*  is  not  strongly  developed . 

Much  changes  when  two  trans  ligands  in  square-planar  ML*  are  bent 

back  to  give  the  angular  fragment  Q.  The  LUMO  becomes  a  hybrid  more 

localized  on  the  metal,  and  resembling  more  the  LUMO  of  ML*  .  And  most 

importantly  the  yz  orbital,  the  one  which  lies  in  both  the  plane  of  deformation 

and  the  plane  of  Ha  approach,  is  destabilized  (thus  moving  it  closer  in  energy 

to  ou*  of  Ha.  with  which  it  interacts)  and  hybridized  away  from  the  fragment 

(thus  providing  better  overlap  with  ou*  of  Ha).  It  is  no  wonder  the  stabilization 

shown  in  Figure  6  is  so  great  -  one  has  moved  part-way  toward  the  product 

geometry  of  oxidative  addition.  The  deformation  toward  an  angular  fragment 

has  activated  ML*  for  the  reaction,  a  point  noted  and  discussed  in  detail 

10 

by  Sevin  and  by  Dedieu  and  Strich.  We  will  soon  see  the  relationship  of 
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this  phenomenon  to  the  observation  of  CH  activation  by  CpML  intermediates. 

We  have  also  calculated  potential  energy  curves  for  the  approach  of 
methane  to  Rh(CO)*'-.  The  perpendicular  approaches  resemble  those  of  Hs  . 

The  parallel  ones  are  dominated  by  steric  repulsion,  so  much  so  that  even 
the  very  attractive  approach  of  Ha  to  angular  ML*  is  transformed  into  a  re¬ 
pulsive  one  in  the  case  of  methane. 

The  observation  of  short  M-H  contacts  in  complexes  of  type  1J3  seems 
in  disagreement  with  our  repulsive  curves  of  Figure  6.  One  could  argue  that 
this  is  the  consequence  of  the  choice  of  carbonyl  ligands  in  our  model  calcula¬ 
tions.  Complexes  of  the  type  always  carry  donor  ligands  such  as  phosphines 
or  halides.  So  we  did  the  same  calculations  replacing  the  four  carbonyls  by 
chlorides.  This  substitution  in  planar  ML*  fragments  is  well  known.  The 
main  result  is  a  moderate  destabilization  of  the  tag  set  (xy,  xz,  yz)  and  a 
large  destabilization  of  z,  which,  however,  becomes  85%  localized  on  the 
metal.  The  z2  orbital  remains  unchanged.  As  the  interaction  between  ML* 
and  Ha  or  CH*  is  dominated  by  the  z2a^2  repulsion,  the  four  E=f(L)  curves 
remain  still  repulsive,  even  if  they  are  less  so. 

A  careful  examination  of  the  structures  of  type  18  show  that  their  M-C 
chains,  because  of  their  steric  encumbrance  and  their  partial  rigidity  are  forced  to 
lie  in  a  plane  roughly  perpendicular  to  the  ML*  plane,  bringing  a  C-H  bond 
in  proximity  to  the  metal.  The  compounds  will  minimize  M***C-H  re- 
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pulsion  by  bringing  the  H  atom  into  an  axial  position  (leading  to  a  positive 

M..  .  H  overlap  population)  and  the  C  atom  as  far  away  as  possible  (thus  minimizing 

the  M***C  negative  overlap  population).  Our  calculations  on  the  model  ^3 


23 

give  overlap  population  values  of  +0.025  for  Rh-H  and  -0.094  for  the  Rh-C, 
the  corresponding  interatomic  distance  being  respectively  1.78  and  2.75  A. 

The  total  M-H-C  overlap  population  is  negative,  in  agreement  with  a  corres- 
r  nding  destabilizing  interaction.  Our  conclusion  is  in  agreement  with  the 
structural  work  of  Echols  and  Dennis'^  who  suggested  that  in  the 
planar  pyrozolyl  Ni  complexes,  there  is  no  M'*’H  short  contact  if  the  mole¬ 
cule  is  sterically  free  to  avoid  it. 

Our  next  goal  is  to  understand  the  activation  of  methane  by  CpML 
intermediates,  such  as  are  thought  to  be  created  in  the  studies  of  Bergman, 
Graham,  Jones  and  their  coworkers.  From  a  theoretical  point  of  view  the 
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irolobal  analogy  d®CpM/L,  24*g“>25  is  obvious,^4  but  does  it  in 

fact  hold  up? 


24  25 

To  answer  this  question  we  constructed  a  hypothetical  reaction  coordinate  for 
the  oxidative  addition  of  CH4  to  Rh(C0)4*  and  CpRhCO,  26.  We  are 


26 


not  able,  as  we  said  above,  to  calculate  a  realistic  path  with  the  extended 
Hilckel  method.  Better  calculations  will  have  to  do  that.  What  we  did  was 
to  make  a  linear  transit  between  a  point  on  the  perpendicular  approach 
(Rh-H  1.755  A,  corresponding  to  the  minimum  on  the  E(L)  curve)  and  an 
idealized  octahedral  product  geometry  (Rh-H  1.6  a,  Rh-C  1.95  A).  The 
transit  is  certainly  not  optimal,  but  does  contain  the  essential  features  of 
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any  reasonable  reaction  coordinate,  for  instance  reorientation  of  the  methyl 
group  to  point  toward  the  metal. 

A  glance  at  Figure  7,  the  comparison  of  the  CpRh  C0(CH4)  and 
Rh(CO)4  +  (CH4)  transits,  shows  how  well  the  isolobal  analogy  works .  Great 
similarity  is  also  seen  in  the  Walsh  diagram  for  the  two  cases.  Given  the 

Figure  7  here 

validity  of  the  analogy  we  would  like  to  go  back  to  the  somewhat  more  sym¬ 
metrical  Rh(CO)4+  model  to  see  how  activation  comes  about. 

A  diagram  showing  the  evolution  of  the  critical  energy  levels  along  the 
reaction  coordinate,  i.e.  a  simplified  Walsh  diagram  for  angular  Rh(CO)4* 
and  CH4  ,  is  shown  in  Figure  8.  The  interacting  fragment  levels  are  not  much 

Figure  8  here 

perturbed  in  the  beginning  of  the  reaction  (left  side  of  Figure  8).  Compare 
with  2£:  at  low  energy  in  ML4  +  CH4  are  xz,  xy  and  z2  ,  higher  is  the 
occupied  yz.  The  final  product  is  a  typical  octahedral  complex  with  an  oc¬ 
cupied  tag  set.  One  of  the  bunch  of  four  highest  occupied  orbitals  is  mostly 
metal  -  H,  CHa  bonding. 

When  we  use  a  square-planar  ML4  fragment  the  calculated  barrier  to 
oxidative  addition  is  substantially  higher.  We  can  trace  the  difference,  i.e. 
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Reaction  coordinate 
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Fig.  7  A  comparison  of  computed  potential  energy  curves  for  hypothetical 
oxidative  addition  reaction  coordinates  of  CH^  to  CpRh  (CO)  and 
Rh  (CO)^.  Both  molecules  are  referred  to  the  same  energy  zero 
when  at  infinite  separation. 
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the  lowering  or  the  activation  energy  when  addition  takes  place  to  the  angular 
fragment  (or  to  CpML)  to  the  higher  energy  of  the  vz  orbital  in  the  latter. 

It  is  worthwhile  at  this  point  to  draw  the  necessary  connection  between 
oxidative  addition  and  its  microscopic  reverse,  reductive  elimination.  The 
latter  reaction  has  been  studied  theoretically  in  some  detail  by  us  and  by 
others.^  The  essence  of  what  happens  in  reductive  elimination  is  that  one 
of  the  two  M-R  z  bond  combinations  goes  down  to  the  new  R-R  0  bond,  27., 
while  the  other  combination  wants  to  correlate  to  R-R  cj*. 


but  cannot,  23.  .  Instead  it  correlates  to  a  metal  orbital.  We  are  describing 
in  words  the  avoided  crossing  so  clearly  visible  at  the  top  of  Figure  8  and 
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reproduced  schematically  in  2_9  below. 


It  now  becomes  clear  that  the  higher  the  metal  vz,  the  higher  the 
heat  of  the  reductive  elimination  (left  to  right  in  29),  but  also  the  lower 
the  activation  energy  for  oxidative  addition  (right  to  left  in  29)  .  The  dif¬ 
ference  between  the  square-planar  and  angular  ML*  fragments  lies  in  the 
energy  of  the  yz  orbitaL  CpML  is  perforce  related  to  the  angular,  and  not 
to  the  square-planar  ML*  fragment,  since  a  cyclopentacienyl  must  elec¬ 
tronically  and  sterically  be  the  equivalent  of  three  facial  and  not  meridional 
ligands . 

For  the  ML*  case  we  calculate  an  activation  energy  for  oxidative  ad¬ 
dition  of  0.92  eV.  This  value  can  be  reduced  farther  if  a  different  starting 
point  for  the  transit  is  chosen. 


* 
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Given  the  importance  of  a  high-lying  yz  orbital  we  can  try  to  think 

of  other  coordination  geometries  that  enforce  such  a  circumstance .  The  d10  ML3 

system  comes  to  mind.  If  it  is  forced  to  be  pyramidal,  as  in  2Q&  or  .£?,  then  it 

12  16 

possesses  a  high-lying  occupied  xz,  yz  pair,  '  propitious  fora  low  ac- 


M 

a  b 
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tivation  energy  for  oxidative  addition.  Our  calculations  on  a  Rh(benzene)- 
+  CH4  model  give  a  barrier  of  0.74  eV,  starting  the  transit  from  an  M-H 
of  1.6  A.  Another  system  with  a  high-lying  yz  level  is  d10ML2,  30c. 

A  calculation  on  Rh(CO)2”  addition  also  gives  a  low  activation  energy.  We 
think  that  such  complexes  merit  investigation  as  possible  C-H  activating 
systems . 


1 
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Metals  and  Metal  Surfaces 

In  principle  we  could  build  up  toward  a  metal  surface  slowly,  by 
examining  a  progression  of  clusters  of  increasing  nuclearity,  31. 
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As  one  does  this  the  levels  multiply  quickly.  We  strain  to  find  frontier  orbital 
arguments  in  which  all  or  most  of  the  responsibility  for  some  basic  act  of 
chemical  reactivity  is  placed  on  one  or  a  subset  of  frontier  molecular  orbitals. 
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The  perturbation  theory  based  language  obviously  remains  valid,  it  just  seems 
that  one  is  doing  more  work  than  necessary  to  trace  down  the  important  in¬ 
teractions,  There  must  be  a  way  of  thinking  about  chemical  reactivity  and 
structure  for  infinite  extended  one-,  two-,  and  three-dimensional  materials 
which  deals  from  the  beginning  in  the  properties  of  bunches  of  levels  and  not 
discrete  levels. 

Such  ways  of  thinking  exist,  of  necessity  cast  in  the  language  of 
17 

solid  state  physics.  That  language  is  not  too  difficult  to  learn,  for  in 
fact  most  of  the  concepts,  though  bearing  different  names,  have  a  one-to- 
one  correspondence  with  constructs  familiar  in  theoretical  chemistry.  So 
instead  of  levels  of  different  point  group  symmetry,  one  has  bands  of  levels 
distinguished  by  a  translational  symmetry  label  which  happens  to  be  a 
vector  in  reciprocal  space,  in  k  space.  The  Fermi  level  is  the  HOMO,  etc. 
The  key  to  thinking  about  groups  of  levels  is  the  density  of  states,  DOS(E), 
the  relative  number  of  states  in  a  given  energy  interval. 

A  typical  density  of  states  curve  for  bulk  Ni,  calculated  by  the  same 
extended  Hlickel  method  as  we  use  for  complexes  is  shown  in  Figure  9.  Note 
the  "d  band",  largely  metal  3d,  between  -8  and  -12  eV.  Above  it  is  a 
broad  s  and  p  band,  the  bottom  of  which  penetrates  substantially  into  the  d 
band.  In  fact  at  the  Fermi  level  the  populations  of  the  various  levels  are 


Total  DOS 


Fig.  9  Oensity  of  states  of  bulk  Ni. 
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as  follows:  d  s  p  Thus  the  s  band  is  one  third  filled. 

The  bulk  Ni  density  of  states  is  characteristic  of  the  rest  of  the 
transition  series.  In  Figure  10  we  show  the  d  band  width  and  the  Fermi 
level  across  the  first  transition  series.  The  Fermi  level  falls  slowly  across 


Figure  10  here 


the  transition  series.  Its  calculated  value  exceeds  the  observed  work  function 
of  the  metal  by  —3  eV,  which  is  a  typical  error  of  the  extended  Hiickel  pro¬ 
cedure  . 

In  our  discussion  of  the  bonding  of  molecules  to  surfaces  we  have 

found  useful  a  measure  of  bonding  common  in  theoretical  chemistry,  but  not 

often  utilized  in  solid  state  physics.  This  is  the  Mulliken  overlap  population 

between  two  atoms.  In  thinking  about  extended  structures  it  is  necessary  to 

think  about  groups  of  levels,  so  the  relevant  quantity  is  the  overlap  population 

weighted  density  of  states,  i.e.  the  relative  number  of  levels  in  a  given 

energy  interval  weighted  by  the  contribution  these  levels  make  to  the  over- 

2 1 

lap  population  for  a  specified  bond.  We  have  called  these  curves  COOP 
curves  (for  Crystal  Orbital  Overlap  Population).  Their  integral  up  to  a  spec¬ 
ified  Fermi  level  is  the  total  overlap  population  for  the  given  electron  count. 

In  preparation  for  our  use  of  COOP  curves  let's  show  one,  for  Ni-Ni 
bonding  in  bulk  Ni,  Figure  11.  The  bottom  of  the  d  band  is  metal-metal 
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Fig.  10  3d  band  width  and  position  of  the  Fermi  level,  as  computed 

It 

by  the  extended  Huckel  method,  for  the  first  transition  series 
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Figure  11  here 

bonding,  the  top  metal-metal  antibonding .  Similarly  for  the  s,  p  band. 

This  is  all  as  expected.  The  total  Ni-Ni  overlap  population  is  0.107  at 
the  Fermi  level. 

We  now  move  to  the  Ni  surface.  The  one  we  have  chosen  is  (111), 
a  closely  packed  surface.  To  take  full  advantage  of  translational  periodicity 
we  must  in  fact  take  a  film  or  slab  of  finite  thickness,  a  typical  tactic  in  solid 
state  theoretical  approaches  to  surfaces.  The  thickness  of  the  film  should  be 
such  as  to  ensure  that  it  approximates  a  real  surface,  yet  also  the  thickness 
must  be  kept  small  for  reasons  of  computational  economy.  Appendix  2  details 
our  studies  of  films  varying  in  thickness  from  a  monolayer  to  5  layers,  and  how 
we  settled  on  3  or  4  layers  as  a  reasonable  approximation  to  a  surface.  We 
have  used  a  four  layer  film  in  our  studies  of  Hs  chemisorption,  and  a  three 
layer  film  for  CH*  activation.  The  four  layer  film  is  shown  schematically  in  32. 
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The  four-layer  film  has  two  identical  surface-like  layers  and  two 
identical  inner  or  bulk-like  layer.  Figure  12.  shows  the  total  density  of  states 
and  its  projection  or  partition  among  the  surface  and  inner  layers.  Note  that 

Figure  12  here 

the  states  in  the  surface  layer  have  somewhat  less  dispersion,  i.e.  form 
narrower  bands.  This  is  true  for  both  the  d  and  the  s,  p  band.  The  reason 
for  this  is  simply  that  the  surface  atoms  have  less  nearest  neighbors  (9)  com¬ 
pared  to  the  inner  atoms  (12).  The  number  of  nearest  neighbors  affects  the 
number  of  interactions  or  overlap  available  to  an  atom,  and  it  is  the  overlap 
which  eventually  controls  the  band  width. 

This  trivial  argument,  summarized  very  schematically  in  33.,  has  a 
nontrivial  consequence.  If  we  plot  for  a  given  band  structure,  say  that  of 
Ni,  the  relative  number  of  electrons  as  one  proceeds  to  fill  the  band,  it 
is  clear  that  the  "bulk"  like  layer  will  fill  first.  The  bulk  will  be  negative 
relative  to  the  surface.  Then  at  about  half -filling  the  two  layers,  bulk  and 
surface  will  be  equally  filled.  As  we  fill  past  this  point  the  surface  will 
become  negative.  The  schematic  plot  of  relative  electron  distribution  as 
a  function  of  electron  count  is  shown  in  34.,  and  is  supported  by  our  com¬ 


putations  . 


Projection  of 
surface  atoms 


Projection  of 
inner  atoms - 


-14 

-16 

-18 

-20 
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X 


Density  of  States  Density  of  States 


Fig.  12  Those  parts  of  the  total  density  of  states  of  a  four-layer  Ni 
(111)  slab  which  are  on  the  surface  and  inner  layers.  The 
dashed  line  is  the  total  density  of  states. 
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Negative  charge — — 


In  reality,  as  one  moves  across  the  transition  series  to  provide  the 
variation  in  electron  count,  the  average  energy  of  a  d  electron,  i.e.  the 
center  of  gravity  of  the  bands  in  3^_,  varies  as  well.  The  general  lines  of 
the  argument  remain  -  at  the  right  side  of  the  transition  series  surfaces  should 
be  negative  relative  to  the  bulk,  at  the  left  side  of  the  transition  series  sur¬ 
faces  should  be  positive.  For  the  specific  case  of  our  four  layer  Ni(lll)  slab 
the  surface  atoms  each  carry  a  charge  of  -0.16,  while  the  inner  atoms  are 


The  above-mentioned  crossover  in  the  behavior  of  the  electron  density 

at  the  surface  as  a  function  of  electron  count  has  been  discussed  by  others  prior 

7  23  7 

to  us.  '  The  reader  is  referred  to  the  excellent  papers  of  Shustorovich 

for  a  general  development  of  the  subject. 

We  are  now  ready  to  bring  H3  onto  the  surface. 


i 
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H5  on  Ni(lll) 

We  cover  both  sides  of  the  four  layer  model  film  of  Mi  with  a  mono- 
layer  of  H3  ,  as  shown  schematically  in  3j5_.  There  is  one  Hz  per  surface 
Ni,  i.e.  a  total  slab  stoichiometry  [NuH*]^.  Each  K2  is  "on  top"  or  per- 
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pendicular" ,  above  a  surface  Ni  atom.  The  H-H  distance  is  frozen  at 
0.74  A,  and  the  nearest  H’^’Ni  distance  is  varied,  L.  We  are  studying  here 
as  closely  analogous  a  situation  as  possible  to  the  perpendicular  approach  of 


Ha  to  a  discrete  MLn  complex. 
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The  total  density  of  states  for  L=3.0,  2.5,  2.0  A  is  shown  in  Figure  13 . 

At  3.0  A  separation  we  would  expect  little  interaction  between  substrate  and 
surface,  and  indeed  what  we  see  is  two  sharp  bands  for  the  monolayers  of  H2  c  and 
c*,  superimposed  on  the  film  band  structure.  Mote,  however,  even  here  the 

Figure  13  here 

slight  destabilization  of  the  c*  band  relative  to  the  free  molecule  value.  As 
the  H2  approaches  the  surface  the  c  orbital  band  retains  its  identity,  persis¬ 
ting  in  its  inefficient  interaction  with  the  surface.  cu* ,  on  the  other  hand,  begins 
to  break  up. 

The  projections  of  the  density  of  states  on  cu*  (i.e.  the  fraction  of 
the  DOS  that  is  cu*)  of  Figure  14  clearly  show  the  strong  interaction  of 
-u*  with  the  s,  p  and  d  bands.  At  1=2.0  A  some  fraction  of  c*  (2%)  has 


Figure  14  here 

been  mixed  into  the  s,  p  band,  especially  around  -3  to  -5  eV.  75%  of 
the  c*  is  pushed  up  to  high  energy,  >  8  eV.  What  is  happening  here  is  a 
typical  2  level  interaction,  36.,  except  that  it  is  now  distributed  over  the 
myriad  o:  levels  of  the  H2C*  and  metal  s,  p  bands. 


The  positions  of  free  Ha  o  and  o*  levels  are  marked  on  the 


There  have  been  a  number  of  theoretical  studies  of  the  interaction  of 

7,25 

Ha  with  transition  metal  surfaces,  modelled  either  by  films  or  by  clusters  '  , 

7 

We  have  mentioned  earlier  the  work  of  Baetzold,  Shustorovich  and  Muetterties  , 

25a 

and  here  refer  further  to  the  work  of  Salem  and  Leforestier  and  of  Lundqvist 
2  S  Id 

and  coworkers  .  The  former  authors  stress  the  importance  of  interactions 
with  both  a  and  o*  of  the  substrate,  much  as  we  will  do  below.  Lundqvist 
has  carried  out  a  thoughtful  analysis  of  metal  surface  adsorbate  interactions. 


In  his  model,  the  H3  bond  breaks  because  the  surface  -  H2  resonance  cor- 
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resonance  corresponding  to  au*  falls  in  energy  and  is  filled  as  H3  approaches 
the  surface.  We  do  not  get  such  an  effect,  but  perhaps  the  analogue  in  the 
extended  Kucha  1  model  is  the  lower  group  of  states  induced  by  c*  in  the  s,  p 
band  as  seen  in  Figure  14. 

We  have  shown  what  happens  to  the  energy  levels.  What  about  the  bond¬ 
ing?  We  can  look  at  the  COOP  curves  as  a  function  of  L  (Figure  15).  At 
large  L  the  H-H  bonding  is  picked  up  as  one  sweeps  through  the  cg  level, 

Figure  15  here 

and  H-H  antibonding  (more  than  the  bonding!)  as  one  passes  through  cu*. 

As  one  moves  the  Ha  in  these  features  persist,  but  now  Ni-H  bonding  enters. 
Look  at  L=2.0  A.  There  is  a  little  Ni-H  bonding  as  one  sweeps  through  the 
Ha^  bond.  Above  that,  until  one  enters  the  main  density  of  cu*,  one  passes 
through  a  region  of  Ni-H  bonding  (and  mirrored  H-H  antibonding .) 

The  sign  of  these  bonding  interactions  is  a  vital  clue  to  the  role  of  o 
and  c*  in  the  bonding  to  the  surface.  There  are  two  extreme  possibilities; 

(a)  c-metal  interaction  dominant;  (b)  c*-  metal  mixing  predominating.  Each 
possibility  has  different  consequences  for  M-H  and  M-M  bonding,  as  we 
will  now  show. 

Suppose  o-M  mixing  were  dominant.  Then  in  some  localized  orbital 
scheme  we  would  get  37.,  simple  in-  and  out-of-phase  mixings  of  c  and 


antibonding  bonding — —  — — antibonding  bonding——— — anlibonding  bonding 
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some  appropriate  symmetry  metal  orbital.  The  resultant  COOP  curve,  when 


this  kind  of  thinking  is  extended  to  a  band  picture ,would  have  the  c  band  K-K 
and  M-H  bonding,  whereas  the  corresponding  metal  band  is  still  H-H  bond¬ 
ing,  but  M-H  antibonding  (38) . 

Suppose  instead  c*-M  mixing  were  dominant.  Then  the  localized  pic¬ 
ture  is  like  39_,  and  the  expected  COOP  curve  4J3.  Now  the  lower,  primarily 
metal  band,  is  M-H  bonding,  but  H-H  antibonding,  through  the  admixture 
of  c*.  The  c*  band  itself  is  both  H-H  and  M-H  antibonding. 


Reality  is  the  superposition  of  the  two  effects.  There  is  no  question 
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what  will  dominate  in  the  region  of  the  bands  derived  from  H3  d~  and  dj*. 
3u:  in  the  region  of  the  bands  derived  from  the  metal  the  two  models  give  op¬ 
posite  predictions:  if  c  mixing  were  dominant  over  c*,  the  intermediate 
region  should  be  H-H  bonding,  M-K  antibonding.  The  reverse  should  be 
true  if  o*  mixing  were  dominant.  Figure  15  gives  a  clear  answer:  In  the 
intermediate  region,  in  the  metal  bands,  the  mixing  of  metal  orbitals  is 
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largely  in  a  bonding  way  with  Ha  ,  and  H-H  bonding  is  weakened  in  the 
same  region.  Clearly  metal  surface-hydrogen  c*  mixing  is  dominant. 

Why  is  this  so?  There  are  both  energy  and  overlap  reasons.  au*  of 
H2  lies  in  the  s,  p  band  to  begin  with.  Second,  the  ou*  coefficients  are 
greater,  so  their  overlap  with  appropriate  symmetry  surface  states  is  perforce 
greater  than  that  of  Cg*. 

So  far  we  have  looked  at  the  most  informative  overall  picture,  independent 
of  electron  count.  But  now  it  is  time  to  focus  our  attention  on  Ni,  and  what 
happens  below  its  Fermi  level,  for  its  particular  electron  count. 

The  overall  charge  flow  and  population  analysis  changes  as  a  function 
of  L  are  given  in  Figure  16.  Ni-H  bonding  is  turned  as  L  decreases,  and 

Figure  16  here 

H-H  bonding  is  decreased.  This  is  accomplished  by  populating  o*  of  Hz  , 
with  relatively  minor  depopulation  of  o.  Note  the  difference  between  ac¬ 
tivation  in  the  discrete  complex  (o*  not  much  populated)  and  on  the  metal 

surface  (o*  reasonably  populated).  The  important  role  cf  the  o*  in  H-H 

„  7,25 

or  C-H  bonding  has  been  stressed  in  a  number  of  previous  studies, 

7 

especially  those  of  Baetzold,  Shustorovich  and  Muetterties. 
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Fig.  16  M-H  and  H-H  overlap  populations  (bottom)  and  a,  a  and  total 
densities  (top)  as  a  function  of  L,  as  approaches  the 
four-layer  slab  of  Ni  (111). 
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To  see  how  the  bond  breaking  occurs  in  derail  we  must  apply  a 
microscope  to  Figures  13-14,  and  zoom  in  cn  the  metal  d  band,  -8  to  -12 
eV.  This  is  where  the  action  takes  place.  Figure  17  shows  the  total  DOS 
of  the  Ni(lll)  four  layer  film  with  Hs  overlayers,  L=2.0  .1,  in  this  smaller 
energy  window.  The  projection  of  these  states  on  z  is  too  small  to  show  up, 
but  that  of  z*  is  clearly  visible.  It  is  the  integral  of  this  projection  up  to 

Figure  17  here 

the  Fermi  level  which  gives  the  0.044  population  of  c*  that  may  be  read 
off  Figure  16. 

So  z*  penetrates  the  d  band  and  is  responsible  for  M-H  bonding. 
But  in  more  detail  how  does  it  do  it?  We  can  apply  a  microscope  to  Figure 
15  and  zoom  in  on  the  COOP  in  the  d  band.  This  is  Figure  18. 

Figure  18  here 

Please  note  that  the  scale  on  the  COOP  curves  is  not  the  same  for  Figures 
18  and  15.  The  Ni-H  and  H-H  curves  mirror  each  othe  r.  Further  insight 
may  be  obtained  by  looking  at  projections  of  the  DOS  of  the  Ni(lll)  film 
alone  on  s,  z3  ,  and  xz ,  yz  components  of  the  surface  layer,  Figure  19. 
These  are  the  prime  orbitals  of  c  and  n  local  pseudo- symmetry  with  respect 

Figure  19  here 


Projected  DOS  of  <r* 


That  part  of  the  total  DOS  (dashed  line)  which  is  the  H2 

(solid  line  hugging  the  energy  axis).  This  is  at  L=2.0 

and  the  energy  window  shows  the  Ni  (111)  d  block  only. 

♦ 

integrated  a  population  is  the  dotted  line. 


>o 


Energy 
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— —  antibondmg  Ponding — — 


Crystal  Orbital  Overlap  Population 


Fig.  18  The  COOP  curve  for  Ni-H  and  H-H  bonding  in  the  d  block  region. 

here  has  approached  to  L=2.0  A.  These  peaks  in  the  COOP 
curve  which  pick  up  maxima  in  the  projected  DOS  curve  of  surface 

p 

s  and  z  are  marked  accordingly. 
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to  the  impacting  Ha  . 

The  features  in  the  CCOP  cur/e  of  Figure  18  clearly  pick  up  corres¬ 
ponding  features  of  the  DOS  of  surface  s  and/or  surface  z2  .  We  have 
marked  the  most  obvious  features  in  the  corresponding  curves.  The  picture 
is  chemically  consistent.  The  surface  interacts  with  the  substrate  Ha  mainly 
through  Ha"..'*  and  surface  layer  s  and  z‘  orbitals. 

So  far  we  have  restricted  ourselves  to  an  "on-top”,  perpendicular  ap¬ 
proach  of  Ha  to  a  N'i  surface.  Clearly  other  sites  of  adsorption  and  the 
parallel  geometry  must  be  considered.  In  face  we  studied  further  five  basic 
geometries  £1  -  £5.  £i_  may  be  described  as  three-fold  perpendicular,  £2 

as  two-fold,  perpendicular,  £3  as  on-top,  parallel,  ££  as  on  a  bond,  para¬ 
llel,  £3.  as  across  a  bond,  parallel.  For  each  of  these  we  studied  a  range  cf 
metal-Ha  separations. 


44 


45 
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It  was  mentioned  earlier  that  we  cannot  trust  the  extended  Huckel  cal¬ 
culations  for  tr.e  energetics  of  bond  forming  or  breaking.  Cne  would  have  liked 
to  be  able  at  least  to  predict  reliably  relative  sites  of  adsorption,  but  unfortunately 
we  cannot  do  that.  All  of  the  approacr.  geometries,  parallel  or  perpendicular, 
give  rise  to  repulsive  energy  curves.  The  method  overestimates  the  four-electron 
repulsive  component  of  the  interaction  energy.  Perhaps  one  could  still  hope  to 
argue  from  the  softness  or  hardness  of  the  repulsion,  but  we  would  prefer  to 
abandon  the  energy  criterion,  and  focus  instead  on  what  we  know  extended  Huckel 
aces  reasonably  well  in  molecules  -  the  method  gives  a  reasonable  estimate  of 
conding  interactions,  especially  those  dependent  on  orbital  symmetry. 

7/e  chose  to  compare  the  various  approaches  of  Ha  at  a  similar  Ni-H 
separation  of  -2.0  .  Table  1  gives  several  calculated  quantities  for  the 

geometries  studied;  the  Fermi  level  energy,  the  change  m  total  energy,  M-H 
and  H-H  overlap  populations  ,  Ha  z  and  populations,  tr.e  chances  in 
populations  of  the  surface  metal  layer,  both  total  and  classified  according  to 
orbital  and  symmetry  type  (s  ,  p0,  pff ,  d<- ,  dff,  d  5) .  and  the  total  change  in 
the  inner  or  bulk-like  layer.  The  changes  are  relative  to  the  free  surface,  and 
in  every  case  the  convention  is  that  a  negative  number  implies  less  of  elec¬ 
tron  density,  or  an  increasing  positive  charge. 


Table  1  here 


Let  us  review  what  happened  in  the  first  on-top,  perpendicular  approach 
The  Fermi  level  is  raised  slightly,  the  approach  is  destabilizing  or  repulsive, 
M-H  forms  and  K-H  begins  to  weaken.  There  is  minimal  effect  on  c,  but 
substantial  electron  transfer  to  of  Ha  .  That  electron  transfer  occurs 
primarily  from  the  surface  layer  and  in  the  surface  primarily  from  the  s  and 
dc(=za)  orbitals.  All  other  effects  are  small. 

Three-fold  perpendicular,  41:  In  this  geometry  there  is  substantial  de- 
stabiliza  cion ,  yet  good  M-H  bond  formation  and  H-H  bond  weakening .  The 
CCC?  curve  of  Figure  20  shows  something  new,  some  features  at  low  energy 
which  indicate  interaction  with  ~  in  addition  to  z*. 

Figure  20  here 

What  is  happening  here  is  that  both  Ha  c  and  c*  can  interact  with 
surface  orbitals  which  are  bonding  between  all  3  metals,  as  shown  schemat¬ 


ically  in  46_. 


hese  orbitals  are  the  3d  orbitals  at  the  bottom  of  the  band,  and 


0 


Crystal  Orbital  Overlap  Population 


Fig.  20  The  COOP  curve  for  Ni-H  and  H-H  bonding  when  Hg 
to  the  surface,  on  top  of  a  three-fold  hollow  si 


is  perpendicular 
te  (geometry  423 
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most  of  the  s  orbitals  in  the  d  band,  since  the  latter  are  the  bottom  of  the 
26 

s  ba  nd . 

The  interaction  is  strong,  because  it  occurs  with  three  metal  atoms  in¬ 
stead  of  one.  a*  still  dominates  the  mixing,  but  there  are  signs  that  o  be¬ 
gins  to  enter  the  picture.  Note  for  instance  the  destabilization,  the  result  of 
(overestimated)  interactions  between  the  surface  and  a.  Also  the  Fermi  level 
rises  slightly. 

One  interesting  feature  which  is  the  result  of  the  stronger  interaction  is 

the  polari2ation  of  the  metal  by  the  adsorbate.  Ha  gains  0.065  electrons, 

but  the  surface  layer  loses  -0.121  electrons.  It  loses  them  to  H2  ,  but  also 

to  the  inner  or  bulk  layer.  What  is  at  work  here  is  interaction  ©  discussed 

in  an  earlier  section,  a  substrate-induced  reorganization  of  electron  density. 

7 


This  interaction  has  also  been  discussed  by  Shustorovich , 
this  process  in  some  detail. 


Let  us  examine 


Suppose  that  there  is  some  distribution  of  levels  in  the  band  such  that 
some  levels  are  more  surface-like  than  bulk-like.  This  is  shown  schematically 
in  47^,  where  the  surface-like  levels  are  marked  as  dashed  lines,  the  bulk- 


47 


i 
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like  levels  as  solid  lines.  If  the  dominant  interaction,  as  tar  as  energy  is 
concerned,  is  repulsive  (in  our  case  with  Ha  a)  then  surface-like  levels 
will  be  pushed  up  mere,  some  above  the  Fermi  level.  The  surface  layer  will 
be  depopulated  relative  to  the  bulk.  And  within  the  surface  those  levels 
primarily  involved  in  interaction  (s,  dc ,  d^)  will  be  depopulated  relative  to 
those  not  participating  in  interaction  (for  example  d*). 

In  the  case  at  hand  we  see  clearly  the  surface  to  bulk  electron  shift, 
and  the  increased  population  of  the  dfc,  in-plane  orbital  set.  This,  in  turn, 
win  cause  some  decrease  in  metal-metal  bonding  in  the  surface,  because  the 
newly  filled  dj-  levels  are  metal-metal  antibonding. 

Two-fold,  perpendicular,  42.:  As  might  have  been  expected,  the  results 
for  this  geometry  are  intermediate  between  the  on-top  and  three-fold  sites. 

QprtoPi  paf^Rel-  43:  The  gross  indicators  of  this  geometry  of  approach 
are  disappointing  -  M-H  is  not  even  slightly  bonding,  and  given  how  good 
a*-  M  interaction  was  in  the  discrete  complex  it  is  startling  to  see  so  little 
population  of  Ka  3*.  How  can  this  be? 

The  solid  state  case  of  course  is  much  more  complicated  than  the  dis¬ 
crete  molecular  one .  For  LaM-Hs  in  a  parallel  geometry  the  Ha  3  orbital 
was  not  allowed,  by  symmetry,  to  interact  with  the  same  tag  orbital  that 
gained  so  much  when  it  mixed  c*  into  itself.  In  the  extended  solid  both  o 
and  c*  of  H3  mix  with  the  metal  t2g  orbitals  everywhere  in  the  interior 
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of  the  Srillouin  zone.  The  surface  analogue  of  the  it  bonding  tsg  orbital  is 
subjected  to  both  stabilizing  and  destabilizing  forces. 

Informative  in  this  respect  is  the  COOP  curve  of  Figure  21  left.  The 
peaks  of  negative  Ni-K  (and  positive  K-H)  correspond  to  peaks  in  the  DCS 
of  s  and  za  .  It  is  clear  that  repulsive  interactions  with  Hs  a  dominate  the 
interaction. 

Figure  21  here 

3ut  there  is  a  glimmer  of  hope.  In  the  lower  region  of  the  band  the 
peaks  of  negative  H-H  overlap  population  (minima  in  the  Ni-H  overlap  pop¬ 
ulation)  correspond  to  peaks  in  the  projected  DOS  of  xz  (Figure  19).  If  the 
■-*  -  xz  overlap  could  be  increased  perhaps  this  interaction  could  be  magnified. 
This  can  be  accomplished  by  bringing  K3  closer  to  the  surface,  .48,  and  by 
stretching  the  Ha  bond,  49 .  Some  fragment  overlaps  which  demonstrate 
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this  are  shown  in  Table  2,  which  also  gives  population  changes  as  a  result 
of  these  distance  changes. 

Table  2  here 

Bringing  H2  closer  to  the  metal  (£8)  increases  the  c*  -xz  overlap 
greatly.  M-H  becomes  bonding ,  and  cr*  significantly  populated .  The  COOF 
curve  shows  increasing  improvement  in  M-Ha  bonding  and  a  growing  role  for 
xz.  Stretching  Ha  makes  for  still  stronger  M-H  bonding  and ,  of  course ,  sub¬ 
stantial  population  of  o*.  The  c*  -  xz  interaction  dominates  the  COOP  curve. 

An  interesting  feature  of  the  COOP  curves  is  that  they  indicate  stronger 
M-Ha  bonding  for  a  total  electron  count  of  9,  i.e.  one  less  than  Ni,  for 
Co.  It  is  this  kind  of  selective  catalytic  information  that  we  hope  to  obtain  in 
future  studies  in  this  area. 

Clearly  we  have  a  dissociative  process  at  work.  But  in  the  early  stages 
of  the  reaction  it  is  repulsive,  and  needs  some  activation  energy.  Another 
parallel  approach  seems  more  promising. 

Qa-.a_--b.Qnd ,  parallel,  44:  There  is  good  bonding  overall,  through  re¬ 
pulsive  contributions  still  dominate .  Each  H3  interacts  significantly  with  two 
metal  atoms.  Characteristic  c  mixing,  shown  in  50 .  is  with  z3  and  ■  in 

the  lower  part  of  the  conduction  band,  and  with  s  throughout  the  band.  There 

is  much  polarization  of  the  metal  and  electron  reorganization  within  the  surface. 
The  Ha  c*  orbital  interacts  predominantly  with  xz  and  z3  combinations  of 
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Table  2.  Interactions  of  Several  Geometries  in  the  On-Top  Parallel  Approach 
of  H3 


Overlap  Populations  Charges 


<  xz  iau*> 

M-H 

H-H 

c  c* 

0 

OrO 

Z  0  A 

0.072 

-0.017 

0.768 

1.985  0.023 

or 

°i 

rO 

• 

1  7d 

0.125 

0.048 

0.732 

1.974  0.042 

on 

»  OA 

O -T-O 

j  1  7  A 

0. 130 

0.186 

0.490 

1.934  0.211 

cb 

i 


-  54  - 


type  5  1 .  in  the  upper  part  c£  the  band.  The  resultant  population  of  c*  is 
good,  and  the  Ni-H  overlap  population  the  largest  in  Table  1.  Once  again 
a  somewhat  lower  electron  countwill  give  better  M-H  bonding. 

The  alternative  across -a-bond ,  parallel  geometry,  45  shows  no 
special  features,  though  its  bonding  interactions  grow  if  H-H  is  stretched. 

To  summarize:  We  cannot  compute  a  potential  energy  surface,  but  with 
some  detective  work  through  the  projected  DOS  and  COOP  curves  and 
populations  we  can  trace  the  origins  of  the  various  oonding  trends.  A  sig¬ 
nificant  aspect  of  this  section  is  our  finding  of  optimum  M-K  bonding  for 
a  parallel  geometry  in  which  Ha  is  lying  over  a  bond,  .44  or  5^.  This 
kind  of  two-metal  assisted  cleavage  of  H2  is  not  possible  in  a  discrete 


i 
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mononuclear  complex. 


H - H 

/ 

m: - m 
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The  important  underlying  theme  of  this  analysis  is  that  the  interaction 
of  Hs  with  a  metal  surface  is  qualitatively  no  different  from  the  similar  in¬ 
teraction  with  the  metal  center  of  a  discrete  molecule.  There  are  important 
differences  in  the  pacing  of  involvement  of  the  Ha  c  and  c*  (c*  is  more  im¬ 
portant  in  the  surface  case).  But  the  fundamental  aspects  of  the  reaction, 
which  we  discussed  in  the  introduction,  remain.  There  is  transfer  of  electrons 
from  Ha  ~ ,  to  H2  c*.  As  a  result  the  H-H  bond  breaks  and  M-H  bonds 
form. 
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Methane  on  Ni(lll) 

A  1:1  coverage  of  a  close  packed  metal  surface  by  methane  is  not 
possible,  no  matter  what  the  approach  geometry.  The  reason  for  this  is 
excessive  steric  hindrance  between  the  methanes.  We  went  to  a  coverage 
of  a  third,  using  the  unit  cell  shown  in  a  solid  line  in  5J3,  instead  of  the 
1:1  coverage  originally  used,  dashed  line.  The  reduced  coverage  is 
(/'3  x/3)R30°  . 

O  O  O  O 


o  o  o  o 
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The  larger  unit  cell  thus  brought  about  forced  us  to  a  three  layer  film  in 

place  of  the  four  layer  one  we  had  used  previously.  We  also  economized 

by  covering. only  one  side  of  the  film. 

Several  geometries  are  feasible,  .54- .58  among  them.  The  methane 

.  27 

molecules,  fixed  tetrahedral  with  C-H~l:l  A  were  placed  so  that  all 

the  closest  Ni-H  contacts  were  2.0  A.  The  various  approaches  are  label- 


i 
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58,  2H  acroas  bond 


led  by  the  number  of  hydrogens  directed  toward  the  surface,  and  by  a  geomet¬ 
rical  descriptor.  Some  calculated  quantities,  paralleling  those  we  found  use¬ 
ful  to  analyze  for  adsorbed  H8  ,  are  given  in  Table  3.  Among  the  three  ta(c) 
orbitals  of  CH4  ,  those  which  are  interacting  preponderantly  with  the  surface 
are  the  ones  bearing  the  largest  coefficients  on  the  interacting  hydrogen(s). 

Table  3  here 
2  6 

These  orbitals  are  shown  in  5.9 . 


Mty  ii.  Ill  I,, yt;,  1  ’V1'1  ‘ 

i •  i  l.lec.uoii 


M1IM  |  •  i  P|1K  •: 1  111  I  »Al  | 


I  H  ,  on  top  IH,  three  fold 


3  H 


2H,  along  bond 


2H,  across  bond 
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Notice  that  in  all  cases,  c  and  c*  have  the  same  phase  relationship  on  the 
interacting  nydrogen(s)  and  consequently  will  interact  with  the  same  surface  or 


bitals.  The  1H  geometries,  on-top  or  over  a  three-fold  site,  are  remarkably 
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similar  to  their  Hs  counterparts.  There  is  repulsion ,  M-H  bonding,  C-H 
bond  weakening,  transfer  of  electrons  from  methane  o  and  into  a*,  loss  of 
electrons  from  surface  s,  and  da,  substantial  polarization  of  the  surface. 

In  general  there  is  more  mixing  with  methane  C-H  a  orbitals.  This  was  ex¬ 
pected,  for  it  will  be  recalled  that  the  CH4  t2  is  ~2  eV  above  Ha  ag ,  and 
thus  closer  in  energy  to  the  metal  d  block.  The  stronger  interaction  with 
methane  o  manifests  itself  in  the  greater  (relative  to  Ha)  depopulation  of 
a,  and  the  greater  reorganization  in  the  metal  (see  the  population  shift  to  the 
inner  layer  and  d&) . 

As  for  H2  ,  5.4,  CH4  in  the  three-fold  site,  gives  more  surface-ab- 
sorbate  bonding  than  .53,  because  of  its  larger  number  of  Ni***H  contacts. 

The  3H  and  2H  geometries,  which  have  several  of  these  contacts, 
are  excellent  for  Ni-H  bond  formation.  In  .56  and  _5_7,  since  the  H  atoms 
are  close  to  on-top  positions,  the  interacting  surface  states  are  s  and  z2  . 

5 8  gives  the  larger  Ni-H  overlap  population,  in  this  geometry  s  surface 
orbitals  are  involved  but  also  interactions  of  type  .60  are  engendered  between 
c*  and  a  piece  of  the  yz  band.  The  geometrical  match  and  overlap  are  ex¬ 
cellent  .  The  culmination  of  double  C-H  bond  breaking  in  this  geometry  would 
be  the  formation  of  surface  hydride  and  methylene. 

Geometries  in  which  one  Ni  atom  is  in  contact  with  2  or  3  H  atoms 
have  been  also  studied.  In  both  cases  negative  Ni  -  CH*  overlap  population 


1 
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are  found,  resulting  from  strong  repulsive  interactions  of  type  G>  These 
occur  via  the  H -atoms  and  also  through  the  C-atom  which  in  these  geom¬ 
etries  is  not  far  from  the  Ni  atom. 

7 

We  have  mentioned  earlier  the  important  recent  study  of  Baetzold 
on  the  interaction  of  hydrocarbons  with  transition  metal  films.  That  work 
is  very  much  in  the  same  spirit  as  ours.  We  are  in  agreement  in  the  direction 
of  charge  transfer  between  methane  and  the  surface,  and  the  importance  of 
the  a*  orbitals.  There  is  some  disagreement  between  our  respective  calcu¬ 
lations  on  the  preference  given  to  different  geometries  of  approach. 

Returning  to  general  considerations,  we  think  that  it  is  interesting  that 
the  overlap  population  corresponding  to  one  interacting  hydrogen  in  5J  is 
about  equal  to  the  Ni-H  overlap  population  in  5^  and  lower  than  the  one 


4 


of  55 .  The  C-H  overlap  population  in  these  three  geometries  is  also  about 
the  same.  This  leads  us  to  think  that  for  cyclic  or  chain  alkanes  chemisorp¬ 
tion  may  proceed  by  a  variety  of  surface  -  H  contacts,  the  major  criterion  of 
stability  being  the  largest  possible  number  of  H  atoms  in  contact  with  the 
surface  that  can  be  produced  by  matching  of  the  surface  and  alkane  geometries. 

Recent  studies  of  cyclohexane  on  Ru(001),  where  contacts  of  kind  >4  are 

, 6f , m  .  ,  .  .  , 

suggested  ,  support  this  idea. 
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Ha  and  CH*  on  Titanium 

The  (001)  surface  of  Ti  h.c.p.  exhibits  hexagonal  packing  very 
similar  to  Ni(lll):  the  stacking  is  of  the  type  ARAB  for  Ti(001)  and 
ABCABC  for  Ni(lll).  The  main  difference  between  the  two  surfaces  can  be 
understood  from  figure  10:  Ni  and  Ti  are  at  each  end  of  a  monotonic  series. 
The  Fermi  level  of  Ti  is  ~2  eV  higher  in  energy  than  the  one  of  Ni.  Con¬ 
sequently  interactions  of  type  ©  metal  acting  as  donor  relative  to  a*  or¬ 
bitals,  are  expected  to  be  greater  with  Ti.  This  is  also  reinforced  by  a  better 
overlap  of  a*  orbital  with  surface  orbitals,  due  to  the  diffuseness  of  Ti 
atomic  orbitals.  On  the  other  hand  the  bottom  of  the  Ti  d  band  is  ~3  eV 
higher  in  energy  than  the  one  of  Ni,  consequently  repulsions  of  type  © 
are  smaller  with  Ti. 

For  H2  dissociation  the  same  geometries  (except  for  4T)  as  for  H2 
on  Ni(lll)  have  been  studied;  the  main  calculated  quantities  are  summarized 
in  Table  4 . 

Table  4  here 

Qualitatively,  there  is  no  big  difference  between  the  bonding  of  H2  on 

Ti(COl)  and  on  Ni(lll).  However,  except  for  Ha  on  top  of  a  metal  atom, 

considerably  larger  metal-H  overlap  populations  are  obtained,  associated 

28 

with  small  H-H  overlap  populations.  Clearly,  Ti  is  much  more  dissoc- 
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iative  than  Ni.  The  geometry  of  type  ^4  is  still  the  best  candidate  for  a 
low  energy  dissociative  process. 

An  approach  configuration  of  type  gives  much  more  metal-H  bonding 
with  the  Ti(001)  surface.  This  is  also  a  consequence  of  the  diffuseness  of  titan¬ 
ium  atomic  orbitals:  The  interaction  of  a*  with  the  nearest  metal  atom  is  of 
fc  type,  this  interaction  is  weak  on  Ni(lll)  but  larger  on  Ti(001),  due  to 
better  overlap  with  diffuse  yz  Ti  atomic  orbitals.  In  addition,  the  inter¬ 
action  with  the  second  nearest  metal  (situated  for  both  surfaces  at  —2,4  -  2.5  A) 
is  negligible  on  Ni(lll)  but  important  on  Ti(001).  The  total  resulting  inter¬ 
action  is  shown  in  projection  in  61.  It  involves  the  bottom  of  the  yz  band. 
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Another  typical  feature  of  Ti(001)  is  that  for  some  geometries. 


namely  .43  and  .44,  the  bulk  layer  is  depopulated.  This  is  the  consequence 


of  another  variant  of  interact 


tion  (aj, 


the  "reverse"  of  47,  shown  in  62, 


62 

Some  empty  surface  states  situated  close  to  the  Fermi  level  are,  by  interaction 
with  c*,  pushed  down  below  the  Fermi  level  and  become  populated,  taking 
their  electrons  from  non-interacting  (bulk  and  surface  d*)  filled  levels  situa¬ 
ted  near  the  Fermi  energy. 

For  methane  on  Ti  the  results  summarized  in  Table  3  are  not  qualitatively 
very  different  from  those  concerning  CH*  on  Ni(lll).  Geometries  .55  and 
58  give  the  strongest  Ni-H  bonds.  As  for  H3  ,  the  metal-hydrogen  overlap 

Table  5  here 

population  is  larger  for  Ti(001)  than  for  Ni(lll)  and  the  H-C  overlap  population 
lower.  One  may  expect  more  dissociative  chemisorption  of  saturated  substrates 
on  the  surface  of  metals  situated  on  the  left  side  of  the  periodic  table. 


i 
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Tables.  Titanium  (001)  -  CH^  Interactions 


Structure 


Overlap 

Populations 

Ti  -  H3"  !  C-Hb 


I 


CH*  Electron 
Densities _ 

o  t  o* 


0.045  0.770  7.983  0.024 

) 


0.115  0.724  7.997  0.124 


0.167  0.767 


7.960  0.089 


0.108  0.767 


7.975  0.057 


0.157  0.748 


t 

7.998  0.169  i 

i 


In  cases  where  there  is  more  than  one  Ti-H  contact,  the  entry  is  for  all  the 
Ti-H  contacts  towards  one  CH4  molecule,  summed. 

Average  of  bonds  pointing  towards  surface. 


i 


VA. 
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Concluding  Comments 

The  problems  that  the  extended  Huckel  method  has  with  bond  distances 
have  not  deterred  us  from  seeking  and  obtaining  an  understanding  of  the  basic 
features  of  H-H  and  C-H  activation  in  discrete  transition  metal  complexes 
and  on  two  transition  metal  surfaces. 

We  have  learned  much  that  is  specific  along  the  way:  why  an  Hs  adds 
sideways  to  a  16  electron  ML5  center, the  role  of  steric  problems  in  CH*  ap¬ 
proaching  a  metal  center,  how  activation  is  achieved  on  daCpML  fragments, 
and  now  it  might  occur  in  d*°ML3  and  ML2  species, how  Ha  interacts  in¬ 
itially  with  a  Ni(lll)  surface,  and  how  that  surface  differes  in  electron  density 
from  a  similar  Ti  surface,  the  apparent  importance  of  a  two-metal  mode  of 
bond  cleavage  on  the  surface. 

But  what  is  most  interesting  about  our  research,  we  believe,  is  die 
demonstration  that  with  proper  tools  it  is  possible  to  illustrate  the  clear 
and  essential  similarity  between  what  happens  in  a  discrete  complex  and  a 
metal  surface.  Indeed,  how  could  anything  very  different  happen,  for  the 
basic  interactions  are  the  same?  In  the  process  of  breaking  an  K-H  or 
C-H  bond  electrons  must  flow  from  a  a  orbital  to  the  metal,  and  from  the 
metal  to  c*.  The  metal-H  bond  forms  at  the  same  time.  To  be  sure  there 


are  differences  in  the  pacing  of  these  electron  transfers.  In  transition  metal 
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complexes  coordinative  unsaturation  is  essential,  and  with  it  the  initial  stages 
of  reaction  are  dominated  by  c-*M  electron  transfer.  But  for  Ni(lll)  the  sur¬ 
face  is  electron  rich,  the  Fermi  level  is  higher  than  for  a  molecule,  and  it  is 
electron  transfer  from  M  -*  c*  that  dominates  the  early  stages  of  the  reaction. 

The  analytical  tools  we  use  in  this  paper  are  a  density  of  states  analysis, 
the  projections  of  that  density  of  states  on  various  atoms  and  orbitals  -  very 
similar  to  a  gross  atomic  populations  for  a  discrete  complex.  We  introduce  an 
immensely  useful  new  indicator  -  the  crystal  orbital  overlap  populations  or 
COOP  curves  (technically  the  overlap  population  weighted  density  of  states). 
This  is  the  solid  state  analogue  of  a  Mulliken  overlap  population  and  allows 
a  limpid  analysis  of  bond  forming  and  breaking  processes.  These  tools,  along 
with  one  we  did  not  use  in  this  paper,  the  extended  structure  analogue  of  Walsh 
diagrams,  give  chemists  a  language  for  understanding  solid  state  structure  and 


reactivity . 
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Appendix  1 .  Extended  Huckel  and  Geometrical  Parameters. 

Molecular  Calculations  .  Extended  Huckel  parameters  for  all  atoms 
used  are  listed  on  Table  6.  Idealized  geometries  were  assumed  and  standard 
band  lengths  and  bond  angles  were  used.  In  the  MLP(n=2-5)  fragments,  all 

Table  5  here 

LML  bond  angles  were  180“  and  90“  .  The  CpRhCo  fragment  was  bent  with 

the  angle  (CC)  (Rh)(Centroid  of  Cp)  equal  to  the  ideal  value  of  1 2 5 . 3 3 ;  the 

Rh  -  centroid  distance  being  1.85  a.  In  the  Rh(Cr  H&  )"  fragment  the  Rh  - 

centroid  distance  was  1.82  A .  All  HCH  angles  were  assumed  to  be  109.47“. 

The  following  standard  bond  distances  were  used:  M  -  CC  =  1.90  A;  C-C  =1.15 

Rh  -  CK3  =1.95  A;  Rh  -  H  =  1.60  A;  Rh  -  ?  =  2 . 3 0  A ;  Rh  -  C  2  =  2.30  A; 

C(Cp)  -  C(Cp)  =  1.43  A;  C  (benz)  -  C(benz)  =  1.41  C-H  =  1.09  i;  H-H  =  0.74 

The  geometries  of  _17  (L  =  CO;  R=R,=H)  and  .23  were  constructed  from  ideal- 

5a ,  d 

ization  of  experimental  structures  using  bond  distances  given  above. 

Bulk  and  Surface  Calculations.  All  the  calculations  were  of  the  tight 
binding  extended  Huckel  type.  The  same  parameters  as  for  molecular  calculations 
(Table  6)  have  been  used  for  C  and  H. 

The  Hj^i's  of  the  transition  metals  from  Ti  to  Mi  have  been  determined 
by  charge  iteration  on  bulk  metals,  assuming  the  charge  cependence  of  metal 
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Table  6.  Extended  Hiickel  Parameters  Used  in  Molecular  Calculations 


Orbital 

Hii  (eV) 

Ci 

r 

b2 

a 

Ci 

Ca3 

Cr 

4s 

-8.66 

1.70 

4p 

-5.24 

1.70 

3d 

-11.20 

4.95 

1.60 

0.4876 

0.7205 

Fe 

4s 

-9.10 

1.90 

4p 

-5.32 

1.90 

3d 

-12.60 

5.35 

2.00 

0.5505 

0.6260 

Rh 

5  s 

-8.09 

2.135 

5p 

-4.57 

2.100 

4d 

-12.50 

4.29 

1.97 

0.5807 

0.5  685 

P 

3s 

-18.60 

1.60 

3p 

-14.00 

1.60 

Cl 

3s 

-30.00 

2.033 

3p 

-15.00 

2.033 

c 

2s 

-21.40 

1.625 

2p 

-11.40 

1.625 

o 

2s 

-32.30 

2.275 

2p 

-14.80 

2.275 

H 

_l_s _ 

-13.60 

1.30 

. 

_ _ 

_ 

Contraction  coefficients  used  in  the  double  £  expansion 
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Hjj/s  given  by  Gray's  eguations  .  The  A,  B  and  C  iteration  parameters 

were  taken  from  Reference  30.  Experimental  h.c.p,  f.c.c.  and  b.c.c. 

31 

structures  were  used  except  for  Mn  for  which  a  b.c.c.  structure  was 
assumed  with  a  lattice  parameter  determined  by  averaging  those  of  Cr  and 
Fe. 

The  extended  Hiickel  parameters  for  Ti  to  Ni  are  listed  on  Table  7. 

Note  that  they  are  substantially  higher  in  energy  than  the  same  parameters  in 
discrete  molecular  complexes. 

Table  7  here 

Calculations  of  Ha  on  Ni  and  Ti  were  made  assuming  a  two  di  mensional 
slab  of  metals,  four  layers  thick,  with  a  lxl  coverage  on  both  sides  of  the 
slab.  Calculations  with  Ha  on  one  side  give  very  similar  results.  Interactions 
between  H3's  have  been  dropped  to  simulate  a  low  coverage.  The  repeating 

32 

unit  cell  contains  four  Ni  atoms  and  two  Hs  molecules.  A  14  k  point  set 

was  used  in  hexagonal  symmetry;  for  lower  symmetries,  special  points  set 

were  generated  by  symmetry  reduction  of  this  hexagonal  set. 

Calculations  of  CH*  on  Ni  and  Ti  were  made  with  a  slab  of  three 

metal  layers  with  a  (/3  x/3)R30°  coverage  of  CH4  on  one  side  of  the 

slab  only.  The  repeating  unit  cell  contains  nine  Ni  atoms  and  one  CH* 

32 

molecule.  A  5  k  point  set  was  used  in  hexagonal  symmetry,  from  which 
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Table  7.  Extended  Huckel  Parameters  Used  in  Metal  Bulk  and  Surface  Calcu¬ 
lations. 


Orbital 

Hu(eV) 

Cl 

Ca 

Ci 

Ca 

Ti  4s 

-6.3 

1.50 

4p 

-3.2 

1.50 

3d 

-5.9 

4.55 

1.40 

0.4206 

0.7839 

V  4s 

-6.7 

1.60 

4p 

-3.4 

1.60 

3d 

-6.7 

4.75 

1.50 

0.4560 

0.7520 

Cr  4s 

-7.3 

1.70 

4p 

-3.6 

1.70 

3d 

-7.9 

4.95 

1.60 

0.4876 

0.7205 

Mn  4s 

-7.5 

1.80 

4p 

-3.8 

1.80 

3d 

-8.7 

5.15 

1.70 

0.5140 

0.6930 

Fe  4s 

-7.6 

1.90 

4p 

-3.8 

1.90 

3d 

-9.2 

5.35 

1.30 

0.5366 

0.6678 

Co  4s 

-7.8 

2.00 

4p 

-3.8 

2.00 

3d 

-9.7 

5.55 

1.90 

0.5550 

0.6678 

Ni  4s 

-7.8 

2.1 

4p 

-3.7 

2.1 

3d 

-9.9 

5.75 

2.00 

0.5683 

0.6292 

Contraction  coefficients  used  in  double  Q  expansion. 
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special  sets  were  obtained  by  symmetry  reduction  for  lower  symmetries. 

In  all  calculations  the  H-H  distance  is  0.74  A  (unless  otherwise 
specified).  The  C-H  distances  were  set  to  1.1  a  if  they  were  interacting 
with  the  surface,  if  not  they  were  set  to  1.09  A.  HCH  angles  were  109.47° 
The  Ni-Ni  and  Ti-Ti  distances  are  taken  from  Reference  31. 
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Appendix  2.  The  Film  or  Slab  Model  for  Surfaces. 

In  order  to  determine  the  best  compromise  between  time  of  computation 
and  accuracy  of  the  model,  a  study  of  the  dependence  of  surface  electronic 
structure  on  slab  thickness  was  undertaken.  In  Table  8  are  listed  several 
computed  quantities  for  Ni(lll)  sLabs  made  up  of  different  numbers  of  layers. 

Table  8  here 

32 

These  were  carried  out  using  a  hexagonal  30  k  point  set  .  Reasonable 
convergence  is  reached  for  a  slab  of  four  layers.  Electron  densities  of  the 
middle  layer  are  not  very  different  from  those  calculated  for  three  dimensional 
bulk  nickel  (using  a  110  k  point  set);  the  major  difference  is  in  the  3d  pop¬ 
ulation,  due  to  the  negative  polarization  of  the  surface  layer  described  in  the 
text. 

The  overlap  population  inside  the  slab  is  close  to  the  one  obtained  for 
bulk  Ni.  That  the  Ni-Ni  overlap  population  is  largest  on  the  surface  may 
be  explained  by  narrowing  of  the  surface  3d  band  as  shown  in  Compared 
to  the  corresponding  inner  slab  states,  the  bonding  surface  states  are  less 
bonding  and  the  antibonding  states  are  less  antibonding.  For  nickel  almost 
all  the  antibonding  d  levels  are  filled.  As  an  antibonding  level  is  in  fact 
more  antibonding  than  the  corresponding  bonding  level  is  bonding,  the  loss 
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of  antibonding  character  dominates  on  the  surface,  causing  the  increased 
Ni-Ni  overlap  population. 

We  also  have  observed  that  for  the  four  layer  slab  good  convergence  is 
reached  for  the  projected  DOS  of  the  surface  layer  and  of  the  middle  layer. 
Comparison  of  Figures  10  and  12  shows  that  the  projected  DOS  of  the  inner 
layer  of  a  four  layer  thick  slab  resembles  the  total  DOS  of  the  three  dimen¬ 
sional  bulk  Ni. 

Calculations  of  various  geometries  of  Ha  on  slabs  of  3  and  4 
layers  show  that,  even  if  computed  electron  densities  and  overlap  populations 
are  slightly  different,  the  general  trend  obtained  for  a  four  layer  slab  is 
conserved  for  a  three  layer  slab. 


f 


Overlap  Populations 


Number  of 
Layers 

Electron  Densities 

Between  Neighbors 

eF 

Layer  No. 

Total 

s 

P 

d 

On  Surface 

inside  slab 
(averaged) 

(eV) 

1 

1 

10.00 

0.66 

0.14 

9.20 

0.171 

-9.04 

2 

1 

10.00 

0.58 

0.19 

9.23 

0.131 

0.109 

-8.71 

l(surface) 

10.10 

0.63 

0.20 

9.27 

3 

2(middle) 

9.79 

0.60 

0.24 

8.95 

0.134 

0.112 

-8.59 

4 

l(surface) 

2  (middle) 

10.16 

9.84 

0.63 

0.61 

0.20 

0.24 

9.33 

9.00 

0.132 

0.110 

-8.56 

1(  surface) 

10.21 

0.63 

0.20 

9.38 

5 

2 

9.89 

0.61 

0.24 

9.04 

0.130 

0.110 

-8.56 

3  (middle) 

9.80 

0.62 

0.23 

8.95 

3  D  Bulk  Ni 


10.00 


0.62  0.24  9.15 


0.107 


-8.47 
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